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We study theoretically the full distribution of transferred charge in a tunnel junction (or quantum 
point contact) coupled to a nanomechanical oscillator, as well as the conditional evolution of the 
oscillator. Even if the oscillator is very weakly coupled to the tunnel junction, it can strongly 
affect the tunneling statistics and lead to a highly non-Gaussian distribution. Conversely, given a 
particular measurement history of the current, the oscillator energy distribution may be localized 
and highly non-thermal. We also discuss non-Gaussian correlations between the oscillator motion 
and tunneling electrons; these show that the tunneling back-action cannot be fully described as an 
effective thermal bath coupled to the oscillator. 

PACS numbers: 



I. INTRODUCTION 

The possibility to observe the quantum mechanics of 
a macroscopic object has sparked significant interest in 
nanoelectromechanical systems (NEMS), which consist 
of a mechanical oscillator coupled to a mesoscopic con- 
ductor. In recent experiments, the oscillator motion has 
been measured with near quantum-limited precision us- 
ing the conductor as a detector—^, and cooling of the 
oscillator by quantum back-action has been observed 4 . 
In these experiments it is the current noise of the con- 
ductor (i.e. the second moment of current fluctuations) 
that is used to measure position fluctuations of the os- 
cillator. The effect of the oscillator on the current noise 
has been well studied theoretically, leading to new un- 
derstanding of back-action and quantum dissipation in 
NEMS^^^ ' 10 ' 11 . However, much more information lies 
in the full probability distribution of transmitted charge 
through the conductor, or the full counting statistics 
(FCS) 12 . In addition to being of theoretical interest, FCS 
is an experimentally accessible quantity and the third 
moment was recently measured in a tunnel junction 13 . 
Still more information may be gained by considering the 
conditional evolution: given a particular measurement 
history for the current, what can we say about the state 
of the oscillator? 

In this paper, we study the full statistics of a tun- 
nel junction (or quantum point contact) coupled to a 
nanomechanical oscillator, a system recently realized in 
experiment 1 '- 2 -. This system is a prime candidate for mea- 
suring FCS in NEMS, since the intrinsic shot noise can 
dominate over other noise sources making it feasible to 
measure the higher moments. Starting from a micro- 
scopic, fully quantum model, we calculate the FCS of 
tunneled charge as well as the conditional evolution of 
the oscillator, and find several surprises that would not 
be apparent in a study of the noise alone. Despite weak 
oscillator-conductor coupling, we find that the oscillator 
can strongly enhance the third and higher moments of 
the FCS, leading to a markedly non-Gaussian distribu- 
tion. This results from long-lived energy fluctuations in 
the high-Q oscillator, which allow correlations between 



the oscillator motion and tunneling electrons to accumu- 
late up to the ringdown time of the oscillator, overwhelm- 
ing the weak coupling strength and dominating the FCS. 
Further, even though the conductor couples linearly to 
the oscillator position, the oscillator state conditioned 
on a particular measurement of current can be highly lo- 
calized in energy. Finally, we discuss non-Gaussian cor- 
relations between the current and back-action force on 
the oscillator that are not captured by treating the tun- 
nel junction as an effective thermal bath. These correla- 
tions arise from the random momentum kicks imparted 
to the oscillator by individual tunneling electrons, which 
cause one half of the back-action to be correlated with 
tunneling. The non-Gaussian correlations lead to signa- 
tures in measurable quantities such as the current noise; 
thus, measuring the current noise could provide a way to 
distinguish quantum back-action from the effects of an 
equilibrium bath. 

Note that FCS were studied previously in a very dif- 
ferent type of NEMS, a charge shuttl e 14 ! 15 . Conditional 
evolution in NEMS were studied using a quantum optics 
approac h 16 : 17 , but these studies miss key features arising 
in our microscopically-derived model. The average cur- 
rent and current noise of the NEMS studied here were 
addressed in Refs^ i 10 ' 11 ; unlike these works, we present 
an exact solution of the master equation and study the 
FCS. 



II. MASTER EQUATION AND ITS SOLUTION 

The Hamiltonian of the coupled system is H = H osc + 
Hi ea ds + Ht, where H osc describes a harmonic oscil- 
lator of mass M and frequency Q including dissipa- 
tion due to an equilibrium thermal bath at temperature 
Tcr^. Electrons in the leads are described by Hi ea d s = 
J2a k £ kC ! ak c a k ™ eVm, where c a k annihilates an electron 
in lead a — L, R, V is the junction bias voltage, and 
the operator rh counts the number of tunneled electrons. 
Ht describes electron tunneling for the experimentally 
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relevant case of weak oscillator-junction coupling 8 ' ^ 10 ' 11 , 
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where A is the lead density of states, ij describes the de- 
pendence of the transmission phase on the oscillator posi- 
tion x, and Y* is the raising operator associated with m, 
e.g. [m, Y"n — yt. \y e focus on an inversion symmetric 
system in which r\ vanishes 8 -. 

We describe the system using a reduced density ma- 
trix p(t) tracking the state of the oscillator and m, the 
number of tunneled electrons. In the weak tunneling 
regime of interest, the off-diagonal (in m) elements of 
p decouple from the diagonal elements. Since our aim 
is to calculate the statistics of m, we need only consider 
the diagonal elements, p(m\t) = (m\p(t)\m). Treating 
Ht perturbatively and making a standard Markov ap- 
proximation which requires eV 3> M7, we derive a mas- 
ter equation for p(m; t)^. Fourier transforming in m, 



p(x;*) = E; 



1 p(m;t), the equation reads 
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Here, J?o describes the coherent dynamics of the oscil- 
lator, and the total damping and diffusion coefficients 
are 7 = 70 + 7i and D = Do + D\. The coefficients 

70 and D = (M^qM}/2) coth (M}/2Tq) are associated 
with the equilibrium bath (ks — 1), while 71 and D\ 
describe back-action damping and diffusion due to the 
tunnel junction. Taking the electronic temperature in 
the leads to be much less than eV—, these are given by 

71 = Ht?/2ttM and D x = MjiTi, where the effective 
temperature due to the tunnel junction is T\ = eV/2^. 
Note that if we average over m (i.e. set x = 0), Eq. |J5J| re- 
duces to the quantum Brownian motion master equation 
for an oscillator coupled to two effective thermal baths 
consisting of the environment and the tunnel junction^. 
Conversely, tracing over the oscillator degrees of freedom 
yields the generating function for the FCS, 



$( X ; t) = tr [p( X ; t)} = E m e^ m P(m; t), 



(3) 



where P(m; t) is the probability that m electrons have 
tunneled in time interval t. Note that the trace of pit") 
over all degrees of freedom is J2 m tr [/K m ; t)] ~ 1- 

The above model has been used to study the average 
current and noise d 10 ' 11 ; here we present its exact solu- 
tion and use it to study FCS and conditional evolution. 
To work with Eq. ([5]) , we first express the reduced density 
matrix in its Wigner representation, 



W(x,p) = ± f _Jdy(x + y\p\x - y) e' 2 ^. 



(4) 



In terms of the Wigner function, Eq. ([2]) may be written 

d t W(x,p;x;t) = (C c i+C q )W, (5) 

where the evolution is described by two Liouvillian oper- 
ators, 
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with the tunneling rate for the oscillator at x = given 
by T = T$eV/2TTh = 2L>iT 2 /fi 2 T 1 2 . The Liouvillian opera- 
tor C c i describes the effectively classical evolution of the 
system: the first line of Eq. (|6a[) corresponds to a classi- 
cal Fokkcr-Planck equation for the oscillator coupled to 
two effective equilibrium baths, the environment and the 
junction; the second line describes tunneling as a classical 
Poisson process characterized by a rate T(t) that depends 
on the instantaneous oscillator position x(t). In contrast, 
C q accounts for quantum corrections to the effectively 
classical evolution. The ^-dependent terms involving 71 
and Di in Eq. (|6b[) describe conditional damping and 
diffusion; these terms represent back-action that is cor- 
related with tunneling. Conditional back-action arises 
because each tunneling electron imparts a random mo- 
mentum kick to the oscillator, implying that the momen- 
tum kicks are correlated in time with tunneling events, 
and shows that the back-action of the tunnel junction 
is not fully described as an effective equilibrium bath. 
This is discussed in detail in Section [V] Note that when 
we add the back-action terms in C q to those in C c i, we 
find that exactly half of the total back-action is condi- 
tional (i.e. includes the factor e lx ). The other half of the 
back-action is uncorrelated with tunneling, and cannot 
be understood in terms of momentum kicks imparted by 
tunneling electrons. We thus have the surprising conclu- 
sion that even during periods where no electrons tunnel, 
there is still back-action diffusion and damping. Heuris- 
tically, even if no electrons tunnel, we nonetheless gain 
information about the oscillator and therefore there must 
be back-action. Finally, the remaining terms in Eq. (|6b|) 
are also quantum in nature, and arise from the difference 
between two tunneling processes involving absorption or 
emission of a phonon of energy Ml. In particular, the 
last term oc rfxd x W does not vanish when we trace over 
the oscillator degrees of freedom, and thus represents a 
quantum correction to the average tunneling rate 8 -^ (cf. 
Eq. (fTTj) below). The same correction is obtained from 
a direct calculation of the tunneling rate using Fermi's 
golden rule. 

Eq. ([5]) may be solved exactly for the physical initial 
conditions of a thermal oscillator state. Such a state is 
Gaussian, and remains Gaussian under Eq. ((5]) for all 
times. We also take m = at time t = 0, since this 



3 



is when we start counting tunneled electrons. Thus, the 
Wigner function may be written in the form 



W(x,p;x;t) 



V p (x-5:)^ + V x (p-p) 
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where we have scaled all quantities by the natural units of 
the zero-point motion of the oscillator, Axo — \J h/2Mfl 
and Ap = y MMl/2. The state is fully characterized 
at all times by its means, x and p, its variances, V x and 
V p , its covariance, V xp , and its normalization, e^. These 
six Gaussian parameters depend on both \ and t, and 
satisfy simple ordinary differential equations which follow 
directly from Eq. (J5|. First, the means satisfy 



d t x(x;t) = ftp 

+ 7i (e <x - 1) 
dtPix'it) = -Six -ip 
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where we have again scaled the position and momentum 
by Axo and Apo, and defined the dimensionless cou- 
pling strength A = AxqTi/tq. The x _ dependence of x 
and p encodes correlations between the oscillator motion 
and m. For example, one can easily show that the irre- 
ducible correlation between x and the nth moment of m 
is ((xm n )) = (— i)"!^f-- Next, the variances and covari- 
ance (also scaled by Ax , Ap ) satisfy 



d t V x (x;t) = 2QV x , i 
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where T = coth (Ml/2T Q ). Again, the x~ 

dependence of these parameters describes correlations 
between x 2 , p 2 or xp and moments of m. Finally, the 
parameter <fi satisfies 



d t <P 



1 
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(10) 

and is directly connected to the FCS as discussed in the 
next section. 

Eqs. (f51- [TU|) have simple analytic solutions in the limit 
of long times, and may be solved numerically for all 



times to arbitrary precision. Before using the equa- 
tions to study the FCS and conditional evolution, we 
emphasize an important difference from previous treat- 
ments of conditional evolution in NEMS: the evolution 
of the variances is conditional, as seen directly from the 
X-dependent terms in Eqs. ©■ This is in stark contrast 
to the standard treatment where the variances evolve in- 
dependently of tunnelin g 16 ! 17 . This is partly due to the 
conditional back-action diffusion in C q discussed above 
(cf. Eq. (|6bp ). which implies that momentum fluctua- 
tions of the oscillator are correlated with fluctuations in 
m and leads to the conditional terms in Eq. (|9b[) . How- 
ever, we also find conditional terms in Eq. (|9a[) that arise 
from the classical part of Eq. ([5]) described by C c i . This is 
because we start with the linear x-dependence of the tun- 
neling amplitude in Eq. ([T]), and it follows that the tun- 
neling rate has both linear and quadratic x-dependence: 



r(x) = T [l + 2Aa; + AV] - 
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(11) 



Standard treatments of conditional evolution neglect the 
quadratic dependence, which in our case is inconsistent 
with the starting Hamiltoniar* 2 ^. We stress that the 
conditional (i.e. x-dependent) and unconditional (i.e. x~ 
independent) terms in Eqs. ([9]) appear at the same order 
in the coupling strength A; there is no a priori reason to 
keep one effect and not the other. The results presented 
below are contingent on the conditional evolution of the 
variances. 



III. FULL COUNTING STATISTICS 

It follows directly from Eq. ([3]) that the generating 
function for the FCS is given by the Gaussian param- 
eter <p via $(x;i) = e^ x;t ). From Eq. (fTU]) we see that 
if the average and variance of the oscillator were simply 
constants, then tunneling electrons would obey Poisson 
statistics with an effective tunneling rate given by (T(x)}, 
obtained from Eq. (fTTj) . However, the oscillator position 
is correlated with tunneling electrons; this correlation en- 
ters Eq. (fTU]) through the x-dependence of x and V x and 
leads to deviations from Poisson statistics. From Eqs. © 
and (|10p we obtain P(m;t), shown at several times in 
Fig. [TJ Even for weak coupling, i.e. A 2 ((x 2 )) <C 1, the 
oscillator can have a dramatic effect on P(m; t) at what 
we call intermediate times, t* < t < I/7, causing it to 
become highly non-Gaussian. Here, 
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(12) 



and T = D/M^f is the net effective temperature of the 
oscillator due to both the tunnel junction and the thermal 
environment. We emphasize that in the relevant limit 
of a high-Q oscillator, the timescale I/7 is much larger 
than 1 /r and thus many electrons have tunneled even for 
intermediate times. 
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FIG. 1: (Color online). Main: P(m;t) at three times. We 
have shifted m by its average {m(t)) and scaled by its stan- 
dard deviation a m (t). Inset: oscillator enhancement of the 
variance (black dashed) and skewness (red solid) of P(m; t) 
versus time. We have taken T =Ti = 10 4 Ml, A = 1(T 3 , r = 
0.2,7o = lCT 5 a These values are based on the experiment 
in Ref. Q], except that we decreased the intrinsic tunneling 
strength tq to operate in the tunneling regime and increased 
the coupling strength A to clearly illustrate our results. For 
these values, the fractional change in the average current due 
to the oscillator is A 2 ((a; 2 )} = 0.02. 



autocorrelation in the integrand is proportional to an en- 
ergy autocorrelation (up to insignificant rapidly oscillat- 
ing terms). This contribution initially scales as (T/Ml) 2 
and decays on the very slow timescale of the oscillator 
ring-down time, 1/7. Thus, long-lived energy fluctua- 
tions in the high-Q oscillator allow its influence to build 
up, eventually overcoming the weak coupling strength 
and dominating the FCS. This enhancement occurs when 
the last term in Eq. (|13ap dominates the first, requiring 
T(X 2 T/hfl) 2 /j ^> 1. This condition can be satisfied even 
when the oscillator contribution to the average current 
e(T(x)) is small, as the ratio r/7 is typically large (e.g. 
r/7 ~ 10 8 in Ref. Q). Further, this same condition en- 
sures f < I/7 from Eq. (|12p. resulting in non-Gaussian 
FCS over a wide range of times. 

If the condition for enhancement is met, then the ef- 
fect is even greater for higher cumulants. For example, 
Eq. (|13b[) contains a term proportional to {{x 2 (t)m 2 (t))) . 
This leads to an oscillator-dependent term in the skew- 
ness similar to that in the variance, with an additional 
factor of to resulting in an extra factor of TA 2 a; 2 and an 
extra time integral. We obtain a three-time x 2 autocor- 
relation which initially scales as (T/M7) 3 and decays on 
the timescale I/7, compensating for the extra factor of 
weak coupling. In general, we find that the maximum 
enhancement for the nth cumulant is roughly 



The significant modification of the FCS is due to 
the seemingly weak dependence of the current on x 2 
(cf. Eq. (fTTj) ). To see this, it is useful to consider the 
first few cumulants of to. From Eqs. §5§ and (|10[) . these 
satisfy [x is again scaled by Axq) 



«m"(t))) ~ (TX 2 Tt/m) n 



(14) 



d t ((m 2 )) = d t {m) + 2T [2 A ((xm)) + A 2 ((x 2 m))] 
d t {{m 3 )) = d t {m) + 3t\2X [((xm)) + {{xm 2 ))] 



(13a) 
(13b) 
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where all of correlations depend on t. The first term in 
each equation corresponds to Poisson statistics, in which 
all cumulants would be equal to (to). Correlations of x 
and x 2 with m emerge naturally in the cumulants due 
to the ^-dependence of the tunneling rate in Eq. (fTTj) . 
because m(t) — J Q dt'T[x(t')). As indicated in Eqs. (fT3"j) . 
this allows oscillator fluctuations to affect the variance 
and skewness of to; since x and x 2 are positively corre- 
lated with to, the cumulants will be increased by these 
correlations. Similar correlations appear in the higher 
moments. These correlations can strongly affect the cu- 
mulants due to the slow decay of energy fluctuations in 
the oscillator, as we now discuss. 

Consider Eq. (|13ap . From the x 2 term in T(x), 
fluctuations in x 2 will lead to fluctuations in to. 
Thus the last term in Eq. (|13a[) leads to a term 
2T 2 \ i f*f*dt 1 dt 2 {{x 2 (t 1 )x 2 (t 2 ))) in the variance. The 
factor A is small due to weak coupling; however, the x 2 



for times t* < i < 1/7. This can be seen directly from 
Eq. (|lip by assuming that fluctuations in to are domi- 
nated by x 2 fluctuations for this range of times. 

Fig. Q] shows that P(m; t) is skewed only for intermedi- 
ate times t* < t < I/7; for short and long times the dis- 
tribution is nearly Gaussian. The enhancement of cumu- 
lants compared to their Poisson values (i.e. ((m n (t))) = 
Tt with no oscillator) is shown in the inset of Fig. [JJ For 
short times, the effects of the weakly coupled oscillator 
have not yet built up and we obtain the Poisson statis- 
tics of the uncoupled tunnel junction. For long times 
t I/7, the contribution to ((to")) from x 2 fluctuations 
simply scales as t (and not as t n ), as t is now much longer 
than the lifetime of a typical oscillator energy fluctuation. 
For long times the oscillator still enhances ((to™)) by a fac- 
tor (T/j) n - 1 (X 2 T/hn) n over the Poisson value Tt, but 
since each cumulant is proportional to t, P(m; t) tends 
to a Gaussian^. To estimate the timescale t* for the 
buildup of enhanced cumulants, note that significant en- 
hancement will occur when the oscillator contribution to 
the variance in Eq. (|13ap is larger than the Poisson contri- 
bution. From Eq. Q3) this requires {T\ 2 Tt/Ml) 2 > Tt, 
which yields Eq. (fT2| for the timescale t*. 

In the range of times where the FCS is strongly in- 
fluenced by the oscillator, P(to; t) is directly related to 
P(x). For a thermal oscillator at temperature T , we have 

P( x ) = \J |^e" mi2 / 4T with x in units of Aa; . Assum- 
ing that fluctuations of x 2 are the dominant source of 
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large m fluctuations, and using Eq. (fTTj) , we obtain 



P(m; t) oc exp 



ATX 2 Tt 



(15) 



for to ^> Ft. This estimate describes the tail of P(m;t) 
very well for times t* <C t <C I/7. 



IV. CONDITIONAL EVOLUTION 



The effects of the oscillator on the FCS are the result 
of correlations between x 2 and m; we can thus gain fur- 
ther insight by studying conditional dynamics. The joint 
distribution P(x, m;t) — (x\p(m;t)\x) is shown in Fig. [2] 
Consistent with the FCS, for short and long times we see 
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FIG. 2: (Color online). Joint distribution P(x,m;t) calcu- 
lated for the same times and parameters as in Fig. [T] We 
have shifted m by its average (m(t)) and scaled x and m by 
their standard deviations, a x (t) and a m (t). 

only small correlations between x and to. P(x, to; t) is 
most striking at times t* < t < 1/7 due to correlations 
between x 2 and to. 

Eqs. (IMT0|) may also be used to find the conditional 
energy distribution, P(E\m;t) — given a particular mea- 
surement history and value of m(t), what is the oscil- 
lator's energy distribution? In Fig. [3] we see that for 
t* <C t <C I/7, the conditional energy distributions are 
highly non-thermal and localized at the energy required 
to produce the given value of to from Eq. (jTTJ) , with width 
given roughly by T . The ability to obtain information 
about the oscillator's energy distribution using a weakly 
coupled detector is somewhat surprising, and is another 
result of long-lived energy fluctuations in the oscillator. 



V. NON-GAUSSIAN CORRECTIONS TO THE 
EFFECTIVE BATH MODEL 

The effects discussed so far are captured by the effec- 
tively classical Liouvillian operator L c i of Eq. Ne- 
glecting the quantum corrections results in an "effective 
bath" model, where the back-action effects of the tun- 
nel junction are treated as arising from a second ther- 
mal bath coupled to the oscillator, and the oscillator is 
treated as a classical variable which sets the instanta- 
neous tunnelling rate. However, the conditional back- 
action damping and diffusion terms in Eq. (|6b[) lead to 
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FIG. 3: (Color online). Conditional distributions P(E\m;t) 
for the same times and parameters as in Fig. [1] The green 
(thin solid) line shows the unconditional distribution (average 
over m). 



non-Gaussian correlations between the junction current 
and back-action force operators I and F that are not cap- 
tured by the effective bath model. These arise because 
even though tunneling is stochastic and imparts random 
momentum kicks to the oscillator, each momentum kick 
occurs at the same time that an electron tunnels. This is 
completely missed in the effective bath model, as it treats 
the junction as a thermal noise source independent of in- 
dividual tunneling events. For example, using Eqs. © 
to calculate ((a; 2 (£)m(i))) in the long time limit, we find 
an enhancement compared to the effective bath model: 



A((x^t)m(t))) 



Q 



(t 



(16) 



This implies the existence of non-Gaussian correlations 
between the the current and back-action force. A di- 
rect quantum calculation of the non-Gaussian correlator 
{(F(ti)F(t2)I(t 3 ))) using Keldysh path integrals follow- 
ing Ref. [l2| leads to the same non-Gaussian correction 
given in Eq. (fTB"|) . 

The non-Gaussian correlations may be understood in 
terms of a simple model of quantum back-action. Wc 
describe the oscillator-independent tunneling current as 
a sequence of (5-functions, I{t) — e^2^L Q 6(t — t n ), 
where the intervals between the t n are exponentially dis- 
tributed. The back-action force of the junction is then 
taken to be Fit) — Y^Lo^n^(t — t„), where £ n is a 
zero-mean random variable describing the impulse im- 
parted to the oscillator by the nth electron. The same 
sequence of times {t n } appears in both lit) and Fit), 
reflecting the fact that back-action arises from the ac- 
tion of individual tunneling electrons. If we then take 
= (fori /to) S mn , our simple model reproduces the 
non-Gaussian correlations obtained from Eqs. (| MTU|) : we 
also obtain the expected back-action diffusion constant 
D-y. From the size of £ n we see that the typical mo- 
mentum kick imparted by a single tunneling electron is 
given by Ap ~ fori/ 'tq, and not by the Fermi momentum. 
This value for Ap is consistent with the Heisenberg uncer- 
tainty principle, since the sensitivity of a position mea- 
surement scales as Aa; ~ tq/t\. We thus have a simple 
picture for the source of the conditional part of quantum 
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back-action: it arises from tunneling electrons imparting 
random momentum kicks of size set by the uncertainty 
principle. Again we stress that this picture only accounts 
for the conditional half of the back-action damping and 
diffusion; the other half is completely uncorrelated with 
tunneling electrons (cf. Eq. (|6bj) and the discussion there- 
after). We also note that one can derive the conditional 
back-action terms in Eq. (|6b[) directly from this simple 
model, from a corresponding classical master equation in 
which each tunneling event is associated with a random 
momentum kick. 

The non-Gaussian correlations discussed above can 
in principle be detected via the finite-frequency current 
noise in the tunnel junction, Si (to). This may be found 
from the time dependence of ((m 2 )) using the MacDonald 
formula—, 

/•OO 

S I {io) = 2e 2 uj \ dt sin (ujt)dt((m 2 (t))). (17) 
Jo 

Note that the frequency-dependent current noise is ob- 
tained from the particle current fluctuations only. In 
the single junction, tunneling is non-resonant and there 
is no place for charge to build up in the system, so 
displacement currents may be safely neglected 2 ^. The 
time derivative of ((w 2 )) is given in Eq. Q13ap . which 
shows that we need the full, time-dependent correlations 
{(x(t)m(t))) and ((a; 2 (i)m(i))) to calculate the current 
noise. These correlations are calculated simply by taking 
the x-derivative of Eqs. ([51-fTU]). The resulting equations 
are readily solved exactly for the physical initial condi- 
tions in which the oscillator is equilibrated with both 
the environment and the tunnel junction, i.e. ((x 2 )) — 

((p 2 )) = ^ m our un it s - The correlation fnnj leads 
to a peak at to = Q in the noise that is very accurately 
captured by the effective bath model. However, the cor- 
relation (^x 2 rnj) leads to peaks at uj = and to = 2£1 
that show signatures of the non-Gaussian correlations. 
This is especially true in the limit 70/71 ~> 7i/T » 1, 
where the non-Gaussian correlations lead to a doubling 
of the current noise peak at to = 0, and completely sup- 
press the peak at to — 2fl as shown in Fig. [U This limit 
requires a back-action damping rate much smaller than 
the intrinsic damping from the environment, and a back- 
action temperature much greater than the temperature 
of the environment. The first of these conditions is nat- 
ural in experiments and the second has been achieved. 
Note that we still require the intrinsic damping of the 
oscillator to be small. 

Measurements of Si(uj) could thus be used to distin- 
guish the tunnel junction's back-action on the oscillator 



from the effects of a Gaussian uncorrelated noise source 
supplied by an equilibrium bath. However, as seen in 
Fig. 21 in the same limit where the signatures are rela- 
tively large, the peak heights themselves are very small 
compared to the frequency-independent shot noise back- 
ground. For this reason, detecting these signatures in the 
current noise would pose a formidable challenge. Further 




FIG. 4: (Color online). Contribution to the current noise 
spectrum near u> = (left) and u = 2fi (right) from to the cor- 
relation ((a; 2 m)). The full calculation including non-Gaussian 
corrections (blue solid) are compared to the results from the 
effective bath model (red dashed). The contributions are nor- 
malized by the frequency-independent shot noise background. 
We have taken T = Ml/2, Ti = 100M, A = 0.01, r = 0.2 
and 70 = 10 -3 fi in order to approach the limit where the 
non-Gaussian signatures are maximized. 



thought will be devoted to more efficient strategies to de- 
tect the non-Gaussian correlations we have identified. 



VI. CONCLUSIONS 

We have studied the statistics of the experimentally 
relevant NEMS of a tunnel junction coupled to a mechan- 
ical oscillator. We have shown that even if the coupling 
is very weak, long-lived energy fluctuations in the oscil- 
lator allow it dominate the FCS. The oscillator-induced 
enhancement of the third moment of the FCS could be 
observed up to measurement times near I/7, well within 
reach of current experiments. We have also shown that 
the effective bath model is not sufficient to fully describe 
the effects of the tunnel junction on the oscillator. Half of 
the back-action is conditional as a result of the random 
momentum kicks imparted to the oscillator by tunnel- 
ing electrons, and this leads to non-Gaussian correlations 
with signatures in the finite frequency current noise. 

This work was supported by NSERC and CIFAR. 
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